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1. Introduction 

>■ 

00 ' Many authors have attempted to model the star formation history of the 

^D ■ universe (eg Cole et al 1994, Pei and Fall 1995) or equivalently the evolution 

riJ \ of the starburst galaxy population (eg Pranceschini et al 1991, 1994, 1997, 

^) • Blain and Longair 1993, Pearson and Rowan-Robinson 1996, Guiderdoni et 

G\ '. al 1998, Dwek et al 1998, Blain et al 1999). Most of these studies explore 

^^ I a small range of models designed to fit the available observations. In the 

r^ • present paper I report the results of a parameterized approach to the star 

^~-*', formation history of the universe, which allows a large category of possi- 

O ' ble histories to be explored and quantified. The parametrized models can 

^_» . be compared with a wide range of source-count and background data at 

jrt ! far infrared and submillimetre wavelengths to narrow down the parameter 

L^ ' space that the star formation history can occupy. 

Throughout this paper I use the most recent Bruzual and Chariot star 
formation scenarios, with a Salpeter IMF with mass cutoffs at 0.1 and 100 
C^ ! Mq, for which the relevant conversion factors are: 

L{Ha)/LQ = 10^-6^*, (1) 
L{m^^m)/LQ = l{f-^^€^^, (2) 
L(28OOl)/L0 = lO^-^^c/)*, (3) 
L(2OOOl)/L0 = lO^-^V*, (4) 
</>* = 42^ z, (5) 

where (p*,<J3z, are the star and metal formation rates in solar masses per 
year per cubic Mpc and e is the fraction of opt-uv light absorbed by dust 
(Rowan- Robinson et al 1997). A Hubble constant of 50 km/ s/ Mpc is used 
throughout, and f^o = 1 (unless otherwise stated). 
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2. Star formation rate derived from ultraviolet data 

I first estimate the star formation history at uv wavelengths, with no cor- 
rection for the eff'ects of dust. For the local, z = 0, star formation rate, I 
have gone back to the Ha data of Gallego et al (1995) in order to remove 
the reddening correction they have applied. I find a mean reddening correc- 
tion of 4.4, to give a value for (/>* of lO^'^-'^^^^-'^^ Mq /yr / Mpc^ . I have also 
estimated the star-formation rate from a sample of 170 bright (B < 13.8) 
nearby (V < 5000 km/s), optically selected galaxies for which there are 
both 60 /im (S(60) > 0.6 Jy) and large-beam 2000 A observations (Buat et 
al 1987, Kinney et al 1993, Deharveng et al 1994, Meurer et al 1995). The 60 
/im and B-band luminosity functions can be used to verify that the sample 
is representative of all galaxies and allow the effective area surveyed to be 
estimated. The 2000 A luminosity density is found to be 10^^ '^'^ W Mpc~^ 
and the corresponding star formation rate is lO~^'^^^^'^^M0yr~^Mpc~^. 
This is consistent with the revised estimate given above from the Gallego 
et al Ha data and with the estimate given by Lilly et al (1996) extrapolated 
from B-band data, W~'^'^'^ MQyr~^ Mpc~^ , but is a factor of 3 lower than 
the value used by Madau et al (1996). 

Madau et al (1996) have given estimates for the star formation rate at 
z = 2 - 4.5 derived from analysis of Lyman drop-out galaxies in the Hubble 
Deep Field. This analysis appeared to demonstrate that the star formation 
drops off steeply at z > 2. Connolly et al (1997) have used photometric 
redshifts for the brighter HDF galaxies to derive the star formation rate at 
z = 0.5 - 2. Their method involves the use of J, H and K data in addition to 
the HST U,B,V,I data to determine the redshifts. Although this certainly 
greatly improves the reliability of the photometric redshifts, it restricts the 
analysis to < 10% of the galaxies detected in the HDF by HST. 

To try to extend the results of Connolly et al to a wider range of redshift, 
and to test the robustness of the conclusions of Madau et al (1996), I have 
used the photometric redshift method of Mobasher et al (1996) to analyze a 
much larger sample of HDF galaxies, namely all those galaxies with Iab ^ 
29.0 from the catalogue of Williams et al (1996) detected in at least the I 
and V bands, a total of 2438 galaxies once duplicates have been eliminated. 
This method gives redshifts accurate to 13 % in (1+z), although there is 
obviously a possibility of aliasing for galaxies detected in only 2 bands 
(22% of the sample). For these 2438 HDF galaxies I have estimated the 
luminosity functions of the galaxies at rest-frame 2800 A in bins of 0.05 
in logio{l + z), for z < 4. The resulting star formation history is shown in 
Fig 1. Agreement with the results of Lilly et al (1996), Treyer et al (1997) 
and Connolly et al (1997) is excellent over the range z = 0.2-2.5. However 
the derived star formation rates at higher redshifts are considerably higher 
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Figure 1. Star formation history derived from ultraviolet data, with no correction for 
effects of dust. Data points are from Gallego et al (1995), (G), Treyer et al (1998) (T), 
Lilly et al (1996) (L), from the sample of bright local galaxies (section 2) at z = 0.01 
(U), and from the analysis of photometric redshifts of HDF galaxies (open circles). The 
models shown are chosen to fit the CDM predictions of Cole et al (1996) (P,Q) = (1.5, 
5.), and to account for the far infrared and submm counts, (P,Q) = (0.5, 5.5). 



than those of Madau et al (1996) and shows no sign of a decline in star 
formation rate beyond z = 1. 



3. Star formation history from infrared or submillimetre data, 
or from ultraviolet data with correction for the effects of dust 

Figure 2 shows the star formation history derived either from ultraviolet or 
HcK data with correction for reddening, or from far infrared or submillimetre 
data under the assumption that most ultraviolet and visible light from star 
forming regions is absorbed by dust (i.e. e ~ 1). 

In the case of the HDF data the extinction is derived using the grain 
properties of Rowan-Robinson (1992) assuming a (conservative) extinction 
corresponding to E(B-V) = 0.1 at the rest-frame wavelength correponding 
to Xohs = 6000 A, which gives correction factors ranging from 1.5 at z =0.2 
to 2.5 at z = 3.8. Most estimates of correction factors to be applied to high 
redshift star-forming galaxies range from 2-7 (Pettini et al 1997, Meurer 
et al 1998, Calzetti 1998, Steidel et al 1999). Note the good agreement of 
my values derived from the HDF with those of Steidel et al (1999) from 
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Figure 2. Star formation history derived from infrared (IRAS: Saunders et al 1990, 
ISO: Rowan-Robinson et al 1997 , F: Flores et al 1997 ), submillimeter (SCUBA, 850 
/im: Hughes et al 1998), and from ultraviolet data with correction for effects of dust, 
from Gallego et al (1995) (Ga), Treyer et al (1998) (T), both corrected by a factor 2 for 
dust extinction, Gronwall (1998) (Gr) and Steidel et al (1999) (S) and from the analysis 
of photometric redshifts of HDF galaxies assuming E(B-V) =0.1 (open circles). 



surveys for Lyman break galaxies. 

The ISO-HDF estimates of Rowan- Robinson et al (1997) for the redshift 
range 0.4-0.7 have been corrected for the fact that one of the 5 galaxies in 
this range is no longer thought to be a starburst (source 3 in Table 1 of 
Rowan- Robinson et al). This has an almost negligible effect on the total 
star formation rate. Also shown are the star formation rates derived by 
Hammer and Flores (1998) from ISO observations of a CFRS field. 

The shape of the observed star formation history is similar to that 
derived from the uv (Fig 1), but the values are a factor 2-3 higher at most 
redshifts, implying e = 2/3-3/4. For comparison Calzetti (1998) gives a 
range for e of 1/2 to 2/3 for starburst galaxies. 



4. Parametrized approach to star formation history 

To study what constraints on the star formation history can be derived 
from source-counts and background intensity measurements at far infrared 
and submm wavelengths, I present here a parameterized approach to the 
problem, investigating a wide range of possible star formation histories. 
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The constraints we have on the star formation rate, 0*(t) are that: (i) 
it is zero for t = 0, (ii) it is finite at t = to, (iii) it increases with z out to 
at least z = 1 ( and from (i) must eventually decrease at high z). 

A simple mathematical form consistent with these constraints is 

Mt)/Mto) = exp Q{l-{t/to)) {t/tof (6) 

where P and Q are parameters ( P > to satisfy (i), Q > to satisfy 
(iii)). I assume that 0*(t) = for z > 10. 

Equation (6) provides a simple but versatile parameterization of the star 
formation history, capable of reproducing most physically realistic, single- 
population scenarios. 

Given an assumed (P,Q) I then determine the 60 /um luminosity func- 
tion, using the IRAS 1.2 Jy sample (Fisher et al 1995). I fit this with the 
form assumed by Saunders et al (1990) 

r;(L) = a(L/L,)i-"e-0-5[Wio{i+L/i.)H' (7) 

and find that the luminosity function parameters can be well approxi- 
mated as follows: 

a = 0.765-0.04 W 

logioL^ = 8.42 + 0.07 W - 2 logwiHo/ WO) 

where W = 0.825 Q - P, and I have assumed fixed values for a = 1.09, 
C* = 0.027 (Ho/lOO)'^. It is not clear that previous studies have correctly 
taken account of the need to change the 60 fim luminosity function as the 
rate of evolution is varied. The study of Guiderdoni et al (1998) explicitly 
violates the known constraints on the 60 /um luminosity function at the 
high luminosity end and as a result the models predict far too many high 
redshift galaxies at a flux-limit of 0.2 Jy at 60 fun, where substantial redshift 
surveys have already taken place. 

To transform this 60 fim luminosity function to other wavelengths I 
assume that the spectral energy distributions of galaxies are a mixture 
of two components, a starburst component and a 'cirrus' component (cf 
Rowan-Robinson and Crawford 1989). I have used the latest predictions 
for infrared seds of these two components by Efstathiou et al (1999). 

There is no explicit allowance for emission from AGN in the models, 
though it is known that at low fir luminosities, about 10 % of galaxies 
contain AGN, with the proportion increasing towards higher luminosities. 
Provided the far infrared component peaking at 60 ;um which I have iden- 
tified as the 'starburst' component is indeed powered by a starburst (cf 
Rowan-Robinson and Crawford 1989, Rigopoulou et al 1996, Lutz et al 
1998), it is irrelevant for estimates of the star formation rate whether the 
galaxy also contains an AGN. I am implicitly assuming that any correlation 
between far infrared luminosity and the optical/uv/X-ray luminosity from 
an accreting black hole is caused by a common feeding mechanism. 
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I can now predict the counts and background intensity at any wave- 
length and by comparing with observed values, constrain loci in the P-Q 
plane. Fig 3 shows a number of such loci for the case of pure luminosity 
evolution in an Einstein de Sitter model (fi = 1). The uncertainties on 
these loci are typically ±0.5 in Q, except for the 60 fiui counts, where the 
uncertainty is much broader (±2 in Q), with rather weak dependence on 
P for P > 0. The locus shown for N60 is in fact the upper end of the 1-a 
range for the number of sources per ster at 60 ;um at 0.25 mJy. Figure 4,5 
shows integral counts at 850 and 60 /Ltm. Figure 6 shows the spectrum of the 
integrated background radiation, compared with selected models. Because 
of the simplistic nature of the assumed optical seds, I have not attempted 
to show optical counts, but it is of interest to try to get roughly the correct 
balance of optical and far infrared/submm background radiation. 

We see that models consistent with the far infrared and submillimeter 
counts can be found, eg (P,Q) = (0.5,5.5), and that such models also are 
consistent with the observed 850 fiui background. However stronger evo- 
lution is required to fit the 140-350 fira background intensity, eg (P,Q) = 
(0.5,6.5), and such models are inconsistent with the observed counts. The 
CDM model, which can be approximated by (P,Q) = (1.5, 5.0), is consistent 
with the 60 /um counts but nothing else. 

An important constraint on the models is that the total mass of stars 
produced in galaxies should be greater than or equal to the mass of stars 
observed, n^ > 0.006 ± 0mi8{Ho/50)-^ (Lanzetta et al 1996), and that 
it should be less than the total density of baryons in the universe, ri* < 
0.05 ± 0.0I(iJo/50)-2 (Walker et al 1991). The models which fit the counts 
( (P,Q) = (0.5, 5.5) ) and 140-350 /xm background ( (P,Q) = (0.5, 6.5) ), 
give fi* = 0.015 and 0.032 respectively. 

We might expect that the cosmological model could have a significant 
effect on the relationship between predicted counts and predicted back- 
ground intensity, since the latter is sensitive to how the volume element 
and look-back time change with redshift. 

To test this we have explored models with A = 0, for which all the 
required formulae are analytic. Figure 7 shows the P-Q diagram for a model 
with Qo = 0.3, with loci similar to those of Fig 3 plotted. In this cosmological 
model there does appear to be a single star formation history which is 
consistent with both source-count and background radiation data, namely 
one with P = 1.5, Q = 8.2. Predicted counts and background spectrum are 
shown for this model as the dotted curve in Figs 4-6. 

Thus within the framework of the types of star formation history con- 
sidered in section 5, there is a clear sensitivity to cosmological parameters, 
with i^o = 0.3 preferred to Qq = 1- This not to say that more complex 
assumptions about the star formation history, for example the inclusion of 
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Figure 3. P-Q diagram, with loci for models fitting 60, 175 and 850 /im counts (solid 
curves) and 140, 240, 350 and 850 /im background intensity (dotted curves). Counts loci 
correspond to N60(0.25Jy) = 2.0, N175(0.1Jy) = 39 (Kawara et al 1998), log N850(4mJy) 
= 900 per sq deg (Hughes et al 1998). Background loci correspond to i^7i/(140) = 12.5, 
!//^(240) = 11, !//^(350) = 6.0 and !/J^(850) = 0.5 nWlm^jsr (Fixsen et al 1998). 



a dust-enshrouded population not represented at low redshift as in Frances- 
chini et al (1997), might not be consistent with Oq = 1. The full dependence 
on cosmological parameters will be explored in later work. 
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Figure 4- Integral source counts at 850 fim. Data are from Hughes et al (1998), Eales 
et al (1999), Small et al (1997). The 3 models shown are, from bottom at faint fluxes, for 
Qo = 1 and (P,Q) = (1.5, 5.0), (0.5, 5.5) and (0.5, 6.5) (solid curves) and for Qo = 0.3, 
(P,Q) = (1.5, 8.2) (dotted curve). 
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Figure 5. Source counts at 60 /im. Data are from Lonsdale et al (1990) (at 0.2-10 Jy), 
Hacking and Houck (1987) (at 50-100 mjy), Gregorich et al (1995) (higher point at 50 
mjy). Models as in Fig 4. 
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Figure 6. Predicted spectrum of integrated background for same models as Fig 4. Data 
from Fixsen et al (1998) (far ir and submm), Pozzetti et al (1998) (opt and uv). 
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Figure 7. P-Q diagram for Q.o = 0.3 model (section 7), with loci for models fitting 60, 
175 and 850 /im counts (solid curves) and 140 and 850 /im background intensity (dotted 
curves) . 
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